Fragility and robustness of the Kardar-Parisi-Zhang universality class 
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We assess the dependence on substrate dimensionality of the asymptotic scaling behavior of a 
whole family of equations that feature the basic symmetries of the Kardar-Parisi-Zhang (KPZ) 
equation. Even for cases in which, as expected from universality arguments, these models display 
KPZ values for the critical exponents and limit distributions, their behavior deviates from KPZ 
scaling for increasing system dimensions. Such a fragility of KPZ universality contradicts naive 
expectations, and questions straightforward application of universality principles for the continuum 
description of experimental systems. Still, we find that the ensuing limit distributions do coincide 
with those of the KPZ class in one and two dimensions, demonstrating the robustness of the latter 
under changes of the critical exponent values. 
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One of the most powerful concepts in contemporary 
Statistical Mechanics is the idea of universality, by which 
microscopically dissimilar systems show the same large 
scale behavior, provided they are controlled by interac- 
tions that share dimensionality, symmetries, and con- 
servation laws. Being rooted in the behavior of equi- 
librium critical systems [l|, universality has more re- 
cently allowed to describe scaling behavior in far-from- 
equilibrium phenomena [^tJIj such as e.g. stock market 
fluctuations [5|, crackling-noise |6|, or random networks 
[7|. In complex systems like these, universality provides 
an enormously simplifying framework, as significant de- 
scriptions can be put forward on the basis of the general 
principles just mentioned. 

Among non-equilibrium systems, particularly conspic- 
uous cases are those displaying generic scale invariance, 
for which scaling behavior occurs throughout parameter 
space [8[ . Examples include self-organized-critical [9| and 
driven-diffusive systems [10[, or surface kinetic roughen- 
ing [ll| . Actually, a paradigmatic model for rou gh i nter- 
faces, the Kardar-Parisi-Zhang (KPZ) equation 



dth = vV'^h + 2 (^'^)^ + Vi^^-t). 



(1) 



is very recently proving itself as a remarkable instance of 
universal behavior. Here, /i(x, t) is a height field above 
substrate position x G M'' at time i, and rj is Gaussian 
white noise with zero mean and variance 2D. Thus, the 
exact asymptotic height distribution function has been 
very recently obtained for d = 1 13|, |lj| : it is given by 
the largest-eigenvalue distribution of large random ma- 
trices in the Gaussian unitary (GUE) (orthogonal, GOE) 
ensemble, the Tracy- Widom (TW) distribution, for glob- 
ally curved ( flat ) interfaces, as proposed in [l5|, see re- 
views in ly, ll7|. Besides elucidating fascinating con- 



nections with probabilistic and exactly solvable systems, 
these and related results are allowing to show that, not 
only are the critical exponent values common to mem- 
bers of this universality class, but also the distribution 
functions and limiting processes are shared by discrete 
models and continuum equations [l8| , and by diverse ex- 
perimental systems, from turbulent liquid crystals [l9[ to 
drying colloidal suspensions |2C)| . 

In view of the success for one-dimensional (ID) sub- 
strates, a natural important step is to assess the behav- 
ior of the KPZ universality class when changing space 
dimension, analogous to e.g. the experimental change 
from 2D to ID behavior for ferromagnetic nanowires, that 
nonetheless occurs within the creeping-domain- wall class 

2l|. Thus, working with discrete models and the con- 
tinuum equation itself, it has been very recently found 

H that, 



for d = 2, indeed universal distributions also 
control fluctuation phenomena in the KPZ class, provid- 
ing higher-dimensional analogs of the TW distributions. 
Again this makes a strong statement on universality, be- 
yond critical exponent values that were already known 
to be shared by the KPZ equation, many discrete mod- 
els |23[, and some experiments 2J], although much less 
than expected [25| . This fact calls for further experimen- 
tal verification, akin to that recently provided [19|, |20[ for 
the one-dimensional case. 

In this work wc report a fragility of the KPZ univer- 
sality class with respect to space dimension. Namely, we 
show that a family of continuum equations that share 
the same symmetries and conservation laws of the KPZ 
equation, and are accurately described by its asymptotic 
distribution function in d = 1, do not share with it even 
the values of the scaling exponents in d = 2. Analogous 
behavior had been found earlier for discrete models of 
(conservative) surface growth [26|, namely, a change of 



the universality class of a given system with d. Here we 
show that it takes place also for non-conserved dynamics, 
and at the level of continuum equations. Note, this is not 
a change in the universality class as a response to changes 
in appropriate system parameters for a fixed dimension, 
as seen e.g. in the context of Barkhausen criticality j27| . 
The lack of universality that we find points to a serious 
difficulty in the identification of the appropriate univer- 
sality class for experimental systems, as it prevents cur- 
sory use of universality arguments to propose theoretical 
descriptions, stressing the need for physically motivated 
models [25|. This fact should be borne in mind, specially 
in view of the timely interest of experimental validation 
for KPZ universality in the 2D case. 
We consider the following equation [28| 

dth^it) = [uk^^ - K.e)h^{t) + ^T[{Vhf]^ + 7?k(i), (2) 

where k = |k|, and i', /C > 0, 7^ is space Fourier trans- 
form, and h]s_{t) and ?/k(0 ^^^ the k-th modes of the 
height and of the noise fields, respectively. Equation ^ 
perturbs the KPZ equation ^ through the introduction 
of the linear term with coefficient v, where < /Lt < 2. 
This family of equations includes celebrated systems, 
such as e.g. the Kuramoto-Sivashinsky (KS) (take /^ — > 2 
and replace k"^ with fc^) and the Michelson-Sivashinsky 
(/i = 1) equations [28|, that combine pattern formation 
at short time and length scales with asymptotic kinetic 
roughening [291 . For specific /x values, Eq. ^ actually 
describes quantitatively surface growth experiments of 
diffusion-limited processes, such as plasma etching |30{ . 



electrochemical [31| , and chemical vapor [32| deposition. 



Equation ^ complies with the standard symmetries of 
the KPZ class. Namely, it corresponds to non-conserved 
dynamics, is isotropic and reflection invariant in the 
substrate coordinates x, invariant under arbitrary shifts 
h -^ h + const., breaks the up-down symmetry h ^-> — /i, 
and satisfies Galilean invariance [33|. Two additional 
features are to be noted, namely, the non-analytic de- 



pendence 3J] of the linear terms on k, and the morpho- 
logical instability j25j . The former induces non-locality 



of the equation when written in real space |35l . |3 



and the latter breaks scale invariance at short time and 
length scales, which is restored back at larg e scales along 
the dynamics, as in the KS system [11|. Indeed, as 
borne out by numerical [28| and dynamic renormaliza- 
tion group [33| results, the asymptotic behavior of Eg. 



11 



([2]) fulfills the Family- Vicsek (FV) scaling ansatz 
implying that the surface structure factor or power spec- 
tral density, S{k,t) = (|/ik(i)|^), scales at long times 
as S{k,t -^ c«) ^ l/fc^"+'', with a well-defined value 
of the roughness exponent. The crossover wave- vector 
value separating white noise from correlated behavior 
also shows FV behavior as fee ~ t^^/^, leading to scaling 
for the global roughness W{t) (root mean square fluctu- 
ation of the surface height) with time a.s W ^ t^ , with 
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FIG. 1: (Color online). Numerical simulations of Eq. ([T]) (•) 
and Eq. ([2]) for ^ = 3/2 (■) and ^ = 7/4 (♦). Surface struc- 
ture factor (upper row) and roughness (lower row) for d = 1 
(left column) and d = 2 (right column). Solid and dashed 
lines represent power-law behaviors as indicated, which are 
discussed in the main text. All the observables have been 
averaged over 10"^ (10^) different noise realizations for d = 1 
(d = 2), starting from a flat initial condition. Error bar are 
smaller than symbol sizes. All units are arbitrary. 



(3 = a/ z. The values of these critical exponents depend 
on /i, and correspond to the KPZ universality class in d 
dimensions, provided ZKPz('i) < /^ < 2. where z-KPz(d) 
is the KPZ value of the dynamic exponent for the corre- 
sponding dimension. Note that, for the morphologically 
unstable condition v > Q that wc consider, the nonlinear- 
ity is dynamically relevant for any value of /i (even if it 
may not control scaling behavior), while for the morpho- 
logically stable situation {v < 0) scaling in Eg. ([ 2| is con- 
trolled by the linear terms for small fi values [28|, |33|, [37 1 . 



In Fig. [T] wc show numerical simulations of Eq. ^ 
for n = 3/2, 7/4 (wc denote these by NL, for non-local) 
and, as a reference, for the KPZ equation itself, for both 
d ~ 1,2, using a pseudospcctral scheme as in [28|, |31j 
and parameters reported in Table [H Both values of fi 
arc larger than or equal to 2:KPz(f) = f/2, thus for 
d = 1 the behavior is well described by KPZ exponents 
/3kpz(1) = 1/3 and aKPz(l) = 1/2. 

Universal behavior here goes beyond exponent values. 
Thus, using the Ansatz [15| 



h^v^t + sgn{\){Ttfx, 



(3) 



we can measure the fluctuations of the interface around 
its mean value Voot, i.e. x = sgn(A)(ft. — Wtx3i)/(rO^- 
Through F we normalize the variance of this stochas- 
tic process to the variance of the TW-GOE distribution 
(that is 0.638) and we are able to compare them. For 
the estimation of Voo and F we followed the procedure 
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TABLE I: Parameters used for the numerical integrations re- 
ported in this work. NL stands for the non-local models, i.e. 
Eq. ^ with fi equal to 3/2 or 7/4. In all cases u — 1. L is the 
size of the ID domain, or the edge of the 2D square, used for 
simulations in Fig. [T] t* is the time used in the computation 
of the quantities reported in Figs. [2] and [3] See more details 
in [331 . All units are arbitrary. 



described in [18| . Values for these constants are reported 
in Table H] (for more details, see [38[). As clearly shown 
in Fig. [21 the random variable x is time- independent and 
distributed according to the GOE TW distribution for 
the three equations. The solid line in Fig.[2]has been cal- 
culated from the solution of the Painleve II differential 
equation [39| , and normalized according to [l5| . 

Further universal behavior is seen to occur, as in other 
instances of the ID KPZ class [18|, for the two-point cor- 
relation function C{x,t) = {h(xo + x,t)h{xo,t)) — (/i)^, 
that scales as C{x,t) ^ {2Tt)'^^^gi{u), where we consider 
u = {Ax/2)/(2TtS^ with gi(u) the covariance of the 
Airyi process [11,01, and A = (2r/A)^/^ for continuum 
equations (for discrete models, A is estimated from the 
local roughness [18[)- As suggested by results in Fig. [21 
even for the ID non-local models considered here, the 
rescaled two-point correlation function collapses perfectly 
onto gi{u), see Fig. [31 

However, when we increase the system dimension 
to d = 2, a remarkable departure from KPZ scaling 
occurs that depends on the relative values of ^ and 
-Zkpz(2) w 1.61. Thus, while Eq. ^ is still well described 
by KPZ exponent values for "large" /i = 7/4 as deduced 
from Fig. [Tl namely, a ~ 0.39 (compare aKPz(2) ~ 0.39 
|23| ) and z ~ 1.61, the "small" /^ = 1/2 system has the 
same exponent values as for d ~ 1. 

This fact bears important consequences on the contin- 
uum modeling of systems, in particular of an experimen- 
tal type, that are presumably in the KPZ universality 
class. Take the ID case as an example. Eq. ([2]) having 
the same symmetries as the KPZ equation, one might 
postulate the latter as a model description for a given 
experiment. But suppose the actual physical interactions 
lead to the occurrence of morphological instabilities (as 
it happens only too often in surface growth experiments 
25l | ) , in such a way that a better description is provided 
by Eq. ([2]) for /.j = 1/2. This will not change the ID scal- 
ing behavior with respect to KPZ universality, even at 
the level of height distributions or correlation functions. 




FIG. 2: (Color onhne). Height distributions for Eq. d (•, 
d = 1; ♦, d = 2) and Eq. © for ^ = 3/2 (a, d = 1; A, d = 2) 
and fi = 7/4 (►, d = 1; T, d = 2). The variable x is defined 
in the text. The solid blue line is the TW (GOE) distribution 
expected for d = 1 [33|. For ID (2D) simulations P{x) is 
estimated from 2048 (1024) independent runs. Inset: zoom of 
main panel, in linear representation. All units are arbitrary. 



However, if one were able to perform an experiment for 
the 2D generalization of the system, a departure from 
KPZ behavior would be obtained, with the conclusion 
that the universality class of the physical system would 
not be KPZ. One might argue that increasing d for a 
fixed fjL makes interactions more non-local in real space 
35|, and that the present fragility of KPZ scaling is only 
superficial |j|. However, this does not circumvent the 
need, for a given physical system, to assess in detail the 
occurrence of e.g. morphological instabilities and/or the 
range of interactions, in order to argue for the correct uni- 
versality class on a safe basis. In any case, this requires 
going beyond symmetry principles to provide the sought- 
for continuum description. Note that, starting out with 
a higher value of // that leads to KPZ scaling both in 
d ~ 1 and 2, such as /i = 7/4, only pushes departure 
from KPZ scaling up to a higher dimension dj/4, such 
that zkpz (1^7/4) > 7/4, which will occur below the upper 
(if finite) critical dimension dc for the KPZ universality 
class, at which ZKPzidc) = 2. 

As a bonus from the above analysis, we obtain, on the 
contrary, an unexpected robustness of the KPZ universal 
height distribution P{x) for d = 2, see Fig. [21 Namely, 
we obtain that the shape of P{x) is robust to changes in 
scaling exponent values, e.g. between those of the KPZ 
equation and those of Eq. ([2]) for /i = 1/2 and d = 2. 
Thus, P{x) turns out to be shared by systems that belong 
to different universality classes, as defined by exponent 
values. We can note that all the equations studied in this 
work share the same so-called Galilean exponent relation. 




FIG. 3: (Color online). ID height-height correlation function 
for Eq. ^ (•) and Eq. O for ^ = 3/2 (■) and ^ = 7/4 (♦). 
Here, u = x^r/2\/{2rt*f^^ while C{x,t = t*) is measured 
from the same surfaces employed to estimate P(x) in Fig. 
[21 The solid blue line provides the covariance of the Airyi 
process [43]. All units are arbitrary. 



a + z = 2, which characterizes, not only the KPZ fixed 



point [ll| , but also the fixed point that governs the low-// 
"phase" of Eq. ^ [28|, [SSf . This scaling relation is quite 
general and occurs in many stochastic systems, as e.g. 
for first-passage percolation and related problems [41[ . 

In conclusion, we have found a fragility of the KPZ uni- 
versality class with respect to space dimension, when per- 
turbed by morphological instabilities combined with non- 
local interactions, within the experimentally substanti- 
ated family of equations, Eq. ^. Note, an important 
perturbation of the KPZ equation by instabilities that 
also respects its space symmetries, is also known to oc- 
cur in the celebrated (noisy) KS system, which is a local 
equation known to lead to KPZ scaling, both in the ID 



111 and 2D [4^ cases. We recall that earlier results have 
also suggested non-universal behavior for the KPZ class 
in d > 1. E.g. for increasing d, details of the noise dis- 
tribution have been reported to become relevant to the 
dynamics |43| . Important issues remain indeed open with 
respect to the dimensional behavior of the KPZ equation 
and class, like the existence and value of an upper crit- 
or even making mathematical sense 



ical dimension 44 



of solutions for d > 1 17 



Although non-equilibrium universality classes are fre- 
uently expected to be more fragile than equilibrium ones 
a highly non-trivial question is to identify, if exis- 
tent, the type of perturbation that is taking place here 
and assess its actual importance [J|. In our case, even 
if the present fragility might be questioned in view of 
the non-local nature of the perturbation, its occurrence 
is not easily circumvented by the symmetry arguments 
that are usually in use for the theoretical description of 
kinetic roughening phenomena. This seems an impor- 



tant caveat^ especially in view of the current quest for 2D 
KPZ scaling behavior in experimental systems. As im- 
plied by our results, in order to assign a universality class 
to a given system, one would need to explore its behav- 
ior under a change in d. However, for many experimental 
systems modifying the space dimension may be hard to 
achieve without significantly altering the basic interac- 
tions that take place. For instance, basic properties of 
fluid flow can drastically change from a quasi-2D Hele- 
Shaw cell to a 3D system, while keeping all additional 
conditions unchanged [45| . This stresses the need for de- 
tailed modeling of the specific peculiarities of the system 
under study, undermining the promise of universality as 
the main toolbox for kinetically rough systems. An anal- 
ogous situation occurs in the context of pattern forma- 
tion, where Goldstone modes associated with the shift 
symmetry h ^ h + const, prevent the existence of a uni- 
versal amplitude equation [46|. In such contexts, model- 
ing has to be done on a system-specific basis. We note 
that in these cases symmetry arguments can be enhanced 
by multiple scales approaches in order to put forward 
general continuum models that successfully describe ex- 
perimental systems 47[. One can ponder SJ] whether 
analogous generalized approaches would be successful in 
the presence of non-localities and noise. 

We have also found that, somewhat unexpectedly, the 
height distributions associated with KPZ universality in 
one and two dimensions are more general than this very 
same universality class, also applying to systems with 
different scaling exponents, but for which the Galilean 
relation a + z ~ 2 holds. This seems reminiscent of 
the occurrence of the same probability distribution for 
global quantities characterizing many equilibrium and 
non-equilibrium strongly-correlated systems [48| . A sim- 
ilar analogy has been recently noticed in the 2D KPZ 
context from a different perspective [49|. On the other 
hand, the TW distribution is known to occur for many 
other disordered systems, like spin glasses [50[. Further 
inquiry into these issues seems a promising route to im- 
prove our understanding of the role of dimensionality in 
the context of the KPZ equation and universality class, 
and for systems with generic scale invariance at large. 
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